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We argue that QCD belongs to a topologically ordered phase similar to many well-known con- 
densed matter systems with a gap such as topological insulators or superconductors. Our arguments 
are based on analysis of the so-called "deformed QCD" which is a weakly coupled gauge theory, 
but nevertheless preserves all crucial elements of strongly interacting QCD, including confinement, 
nontrivial 9 dependence, degeneracy of the topological sectors, etc. Specifically, we construct the 
so-called topological "BF" action which reproduces the well known infrared features of the theory 
such as non-dispersive contribution to the topological susceptibility which can not be associated 
with any propagating degrees of freedom. Furthermore, we interpret the well known resolution of 
the celebrated 1/(1) a problem when would be r\ Goldstone boson generates its mass as a result 
of mixing of the Goldstone field with a topological auxiliary field characterizing the system. We 
identify the non-propagating auxiliary topological field in BF formulation in deformed QCD with 
the Veneziano ghost (which plays the crucial role in resolution of the U(1)a problem). Finally, we 
elaborate on relation between "string -net" condensation in topologically ordered condensed matter 
systems and long range coherent configurations, the "skeletons", studied in QCD lattice simulations. 



I. INTRODUCTION AND MOTIVATION 

The main motivation for the present studies is some 
recent lattice QCD results which can not be easily inter- 
preted in terms of the conventional quantum field theory 
with a gap. To be more specific, the gauge configura- 
tions studied in [1-4] display a laminar structure in the 
vacuum consisting of extended, thin, coherent, locally 
low-dimensional sheets of topological charge embedded 
in 4d space, with opposite sign sheets interleaved. A 
similar structure has been also observed in QCD by dif- 
ferent groups [5-10] and also in two dimensional CP^ -1 
model [11]. A correlation length of the percolating ob- 
jects is order of size of the system ~ L while the width of 
these objects apparently vanish in the continuum limit. 
This is in drastic contrast with conventional expectation 
that in a gapped QCD in hadronic phase the fluctuations 
should have a typical scale ~ Aq CD while the gauge in- 
variant correlation functions must display a conventional 
exponentially weak sensitivity to the size of the system 
~ exp(— AqcdL). Furthermore, the studies of localiza- 
tion properties of Dirac eigenmodes have also shown ev- 
idence for the derealization of low-lying modes on effec- 
tively low-dimensional surfaces. 

It is important to emphasize that the observed long 
range configurations can not be identified with any prop- 
agating degrees of freedom such as physical gluons with 
transverse polarizations. In particular, these long range 
objects contribute to the contact term in topological sus- 
ceptibility which has an on opposite sign in comparison 
with conventional contributions related to propagating 
physical states, see details below. In different words, lat- 
tice studies are consistent with non-dispersive nature of 
these long range objects. 

It is not a goal of the present paper to cover this sub- 
ject [1-11] with a number of subtle points which accom- 
pany it. Instead, we attempt to understand the observed 
long range structure on a deeper theoretical level using 
"deformed QCD" as a toy model [12, 13]. This is a sim- 



plified version of QCD which, on one hand, is a weakly 
coupled gauge theory wherein computations can be per- 
formed in theoretically controllable manner. On other 
hand, the corresponding deformation preserves all the 
relevant elements of strongly coupled QCD such as con- 
finement, degeneracy of topological sectors, nontrivial 9 
dependence, and many other important aspects which 
allow us to test some fascinating features of strongly in- 
teracting QCD, including the long range order observed 
in lattice QCD simulations. It is claimed [12, 13] that 
there is no phase transition in passage from weakly cou- 
pled "deformed QCD" to strongly coupled real QCD. 
The ground state in "deformed QCD" is saturated by the 
fractionally charged weakly interacting pseudo-particles 
(monopoles) which live in 3 dimensions. 

Precisely 3d feature of this model offers a new perspec- 
tive on possible deep relation between "deformed QCD" 
and previously studied (2 + 1) condensed matter (CM) 
systems with a gap which are known to lie in topolog- 
ically ordered phases [14]. These CM systems include, 
but not limited to such systems as: quantum hall states 
as described in [14-16], superconductors as presented in 
[17], and topological insulators as treated in [18], see re- 
cent review paper [19] with large number of original ref- 
erences on analysis of the topological order in condensed 
matter systems. The crucial element which is common 
for all these CM systems is that they can be formally de- 
scribed using a pure topological Chern-Simons effective 
Lagrangian written in the so-called "BF" form. 

One of the main objectives of this work is to derive 
a similar Chern-Simons theory for the "deformed QCD" 
to argue that it also lies in topologically ordered phase. 
Of course, there is a fundamental difference between CM 
systems defined in Minkowski space and Euclidean 3d 
"deformed QCD". In particular, instead of propagat- 
ing quasiparticles in CM systems we have pseudoparticlcs 
which saturate the path integral. Furthermore, the "de- 
generacy" in "deformed QCD" is related to degenerate 
winding states \n) which are connected to each other by 
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large gauge transformation, and therefore must be iden- 
tified as they correspond to the same physical state. It 
is very different from conventional degeneracy in topo- 
logically ordered CM systems when distinct degenerate 
states are present in the system. Nevertheless, the sys- 
tems in both cases demonstrate a huge sensitivity to the 
large distance physics formulated in terms of the bound- 
ary conditions, see few additional comments in [20] on 
analogies between topologically ordered CM systems and 
long range features in QCD. In what follows we dis- 
cuss some manifestations of this long range order in QCD 
which may be considered as a complementary tool, along 
with well-established characteristics such as topological 
entanglement entropy [21, 22]. 

We start our presentation in section II by reviewing 
the relevant parts of the model [12, 13]. In section III 
we argue that the infrared description can be formu- 
lated in terms of Chern-Simons effective Lagrangian in 
BF form. We reproduce the expression for the topologi- 
cal susceptibility using this BF-type action. In section IV 
we present the physical interpretation of the obtained re- 
sults by using analogy with CM systems. In section V we 
consider some profound consequences of our findings for 
high energy particle physics if topological order persists 
in a passage from weakly coupled "deformed QCD" to 
strongly coupled physical QCD. Section VI is our Con- 
clusion where we propose to test very fascinating topo- 
logical properties of the QCD quantum vacuum in a CM 
laboratory. 



II. DEFORMED QCD 

Here we overview the "center-stablized" deformed 
Yang-Mills developed in [12, 13]. In the deformed the- 
ory an extra term is put into the Lagrangian in order to 
prevent the center symmetry breaking that characterizes 
the QCD phase transition between "confined" hadronic 
matter and "deconfined" quark-gluon plasma. The na- 
ture of the gap in this model is reviewed in section II A, 
while in section II B we review the computation of the 
non-dispersive contact term in topological susceptibility 
[23] . This term will be our starting point in construction 
of the corresponding Chern Simons Lagrangian in section 
III. 



A. The model 

We start with pure Yang-Mills (gluodynamics) with 
gauge group SU(N) on the manifold R 3 x S 1 with the 
standard action 

S YM = I d'x -Ltr [**(*)] , (1) 
and add to it a deformation action, 

AS= f d 3 x-Lp[n(x)}, (2) 



built out of the Wilson loop (Polyakov loop) wrapping 
the compact dimension 



O(x) 



(3) 



Parameter L here is the length of the compactified dimen- 
sion which is assumed to be small. The coefficients of the 
polynomial P [O(x)] can be suitably chosen such that the 
deformation potential (2) forces unbroken symmetry at 
any compactification scales. At small compactification L 
the gauge coupling is small so that the scmiclassical com- 
putations are under complete theoretical control [12, 13]. 

As described in [12, 13], the proper infrared descrip- 
tion of the theory is a dilute gas of N types of monopolcs, 
characterized by their magnetic charges, which are pro- 
portional to the simple roots and affine root a a £ A a ff 
of the Lie algebra for the gauge group U(1) N . For a fun- 
damental monopole with magnetic charge a a € A a g, the 
topological charge is given by 
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and the Yang-Mills action is given by 
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The ^-parameter in the Yang-Mills action can be included 
in conventional way, 
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with F^ u = t WF, 



per- 



The system of interacting monopolcs, including 8 pa- 
rameter, can be represented in the dual sine- Gordon form 
as follows [12, 13, 23], 
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where £ is magnetic monopole fugacity which can be 
explicitly computed in this model using conventional 
scmiclassical approximation. The 9 parameter enters 
the effective Lagrangian (7) as 6/N which is the direct 
consequence of the fractional topological charges of the 
monopolcs (4). Nevertheless, the theory is still 27r peri- 
odic. This 2ir periodicity of the theory is restored not 
due to the 2tt periodicity of Lagrangian (7). Rather, it is 
restored as a result of summation over all branches of the 
theory when the levels cross at 6 = ir(mod 2ir) and one 
branch replaces another and becomes the lowest energy 
state as discussed in [23]. 

Finally, the dimensional parameter which governs the 
dynamics of the problem is the Debye correlation length 
of the monopole's gas, 



ml = LC 



(8) 
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The average number of monopoles in a "Debye volume" 
is given by 



9 \ 3 1 



»1, 



(9) 



The last inequality holds since the monopole fugacity is 
exponentially suppressed, £ ~ e~ x l 9 , and in fact we can 
view (9) as a constraint on the validity of the approxi- 
mation where semiclassical approximation is justified. 



B. Topological susceptibility 

The topological susceptibility x which plays a crucial 
role in resolution of the U(1)a problem [24-29] and is 
defined as follows 1 



X(6 = 0) 



do 2 



(10) 
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where 9 is the 9 parameter which enters the Lagrangian 
(6) along with topological density operator q(x) and 
E va ,c(9) is the vacuum energy density determined by (7). 

It is important that the topological susceptibility x 
does not vanish in spite of the fact that q = d^K^ 1 is 
total divergence. Furthermore, any physical state gives 
a negative contribution to this diagonal correlation func- 
tion 



Xdispcrsivc ~ lim / d A xe ikx (T{q(x),q(0)}) 



(ii) 
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where m n is the mass of a physical state, k — > is its mo- 
mentum, and (0|g|n) = c n is its coupling to topological 
density operator q(x). At the same time the resolution of 
the U(1)a problem requires a positive sign for the topo- 
logical susceptibility (12), see the original reference [26] 
for a thorough discussion, 

Xnon-dispcrsivc = Mm / d 4 xe lkx (T{q(x) , q(0) }) > 0.(12) 

Therefore, there must be a contact contribution to Xi 
which is not related to any propagating physical degrees 
of freedom, and it must have the "wrong" sign. The 
"wrong" sign in this paper implies a sign which is op- 
posite to any contributions related to the physical prop- 
agating degrees of freedom (11). In the framework [24] 
the contact term with "wrong" sign has been simply pos- 
tulated, while in refs.[25, 26] the Veneziano ghost had 



1 We use the Euclidean notations where path integral computa- 
tions are normally performed. 



been introduced into the theory to saturate the required 
property (12). Furthermore, as we discuss below the con- 
tact term has in fact the structure x ~ / d A x5 A (x). The 
significance of this structure is that the gauge variant 
correlation function in momentum space 



lim 

fc-s-0 
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develops a topologically protected "unphysical" pole. 
Furthermore, the residue of this pole has a "wrong sign" , 
which precisely corresponds to the Veneziano ghost con- 
tribution saturating the non-dispersive term in gauge in- 
variant correlation function (12), 

(q(x)q(O)) ~ (d»K»{x), 8^(0)) ~ 8\x) (14) 

The singular behaviour (q(x)q(0)) with "wrong sign" has 
been well confirmed by the lattice simulations [3, 5-7, 30]. 

The topological susceptibility in "deformed QCD" 
model can be explicitly computed as this model is 
a weakly coupled gauge theory. In this model it is 
saturated by fractionally charged weakly interacting 
monopoles, and it is given by [23] 



Xym 



Jd i x{ q (x), q (0)) = ^ I Jd 3 x[S( X )} 



(15) 



It has the required "wrong sign" as this this contribu- 
tion is not related to any physical propagating degrees 
of freedom, and it has a <5(x) function structure which 
implies the presence of the pole (13). However, there are 
no any physical massless states in the system, and com- 
putations [23] leading to (15) are accomplished without 
any ghosts or any other unphysical degrees of freedom. 
Instead, this term is described in terms of the tunnelling 
events between different (but physically equivalent) topo- 
logical sectors in the system. 

One should emphasize that the 6{x) function which 
appears in the expression for topological susceptibility 
(15) is not an artifact of small size monopole- approx- 
imation used in [23]. Instead, this singular behaviour 
is a generic feature which is shared by many other mod- 
els, including exactly solvable two dimensional Schwinger 
model and QCD when the contact term is saturated by 
the Veneziano ghost. In fact, this singular behaviour is 
measured in the QCD lattice simulations at strong cou- 
pling [3, 5-7, 30] as we already mentioned. One should 
emphasize again that the significance of this structure is 
that the gauge variant correlation function in momen- 
tum space develops a topologically protected "unphysi- 
cal" pole similar to eq. (13) which will play an important 
role in our discussions in section III when this pole will 
be identified with unphysical topological fields in BF- 
formulation of the theory. 

The 6{x) function in (15) should be understood as total 
divergence related to the infrared (IR) physics, rather 
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than to ultraviolet (UV) behaviour as explained in [23] 



III. BF THEORY FOR DEFORMED QCD 



X YM ~ I 5(x)d 3 X 

( X ^ L \ 
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In different words, the non-dispersive contact term with 
"wrong" sign (15) is highly sensitive to the boundary con- 
ditions and behaviour of the system at arbitrary large dis- 
tances. Therefore, it is natural to expect that a variation 
of the boundary conditions would change the topologi- 
cal susceptibility (15) despite of the fact that the system 
has a gap (8). We'll reproduce the <J(x) function in (15) 
in terms of the topological quantum field theory for de- 
formed QCD constructed in the next section IIIB. This 
will further illuminate the IR nature of the contact term. 

The light quarks can be easily inserted into the system. 
The corresponding generalization of eq. (15) reads [23] 



Xqcd 



d A x(q(x)q(0)) 
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The first term in this equation has non-dispersive na- 
ture and has the positive sign. As explained above this 
contact term is not related to any physical propagating 
degrees of freedom. Instead, it emerges as a result of 
the tunnelling transitions between the degenerate topo- 
logical sectors 2 . The positive sign of this term is the 
crucial element for the resolution of the U(1)a problem. 
The second term in eq. (17) is a standard dispersive con- 
tribution, can be restored through the absorptive part 
using conventional dispersion relations, and has a nega- 
tive sign in accordance with general principles (11). This 
conventional physical contribution is saturated in this 
model by the lightest rj field. It enters xqcd precisely 
in such a way that the Ward Identity (WI) expressed as 
XQCD{jn q = 0) = is automatically satisfied as a result 
of cancellation between the two terms. If the contact non- 
dispersive term with "wrong sign" was not present in the 
system, the WI could not be satisfied as physical states 
always contribute with negative sign in eqs. (11,17). 



2 In the context of this paper the "degeneracy" implies there exis- 
tence of winding states \n) constructed as follows: T\n) = |n+l). 
In this formula the operator T is the large gauge transformation 
operator which commutes with the Hamiltonian [T, H] = 0. The 
physical vacuum state is unique and constructed as a superposi- 
tion of \n) states. In QFT approach the presence of n different 
sectors in the system is reflected by summation over n £ Z in 
definition of the path integral in QCD. It should not be con- 
fused with conventional term "degeneracy" when two or more 
physically distinct states are present in the system. 



In many respects the deformed QCD which is an effec- 
tive 3d theory defined by the action (7) is very similar 
to the abelian Higgs model describing superconductivity. 
The Higgs theory is known to belong to the topologically 
ordered phase [17]. Therefore, our goal here is to provide 
some arguments suggesting that these two models are in 
fact behave very similarly in the infrared. This similarity 
suggests that the deformed QCD also lies in a topologi- 
cal phase. Separately, it has been claimed [12, 13] that 
the passage from the deformed QCD to strongly coupled 
QCD is smooth, without any phase transitions. If true, it 
would imply that QCD also lies in a topologically ordered 
phase. 



BF theory for the deformed QCD. 
Construction. 



The action (7) describes the theory with a gap, simi- 
lar to the conventional Landau-Ginsburg effective action 
describing the abelian Higgs model. In both cases the 
crucial part related to the topological order is missing 
in these effective actions. We want to reconstruct this 
missing topological term. We follow [17] in this deriva- 
tion. As our 3d model (7) has the Euclidean metric, we 
proceed with the Euclidean path integral computations. 
It is different from physical case of the Minkowski 2+1 
BF theory discussed in [17] when one can discuss the 
braiding phases of quasiparticlcs, the Hamiltonian for- 
mulation, etc. Nevertheless, the crucial points can be 
explained using the Euclidean path integral approach. 

We wish to derive the topological action for deformed 
QCD by using the same technique exploited in ref.[17] 
for the Higgs model. We start with construction of the 
source term for a configuration of AfW point-like mag- 
netic monopoles coupled to <t field entering the low en- 
ergy action (7) 



Sj =-P I d 3 xJ(x)a 

Jm 3 



(18) 



N M (a) 



y 0=1 k=i 



'6(xr> - x)a a , 



where = ±1, and a a € A a ff is the affine root of 

the Lie algebra for the gauge group U(1) N and beta is 
defined as 



L \2tt) 



(19) 



It has been demonstrated in [12, 13] that summation over 
all possible positions and orientations of the monopoles 
leads to the low energy action (7), see also [23] with some 
technical details. 

The er field plays the role of the scalar magnetic poten- 
tial as one can see from expression for the U(1) N mag- 
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netic field B l = e ijk4 F jk /2g 
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Essentially, the a^ a ' fields are the physical photons in 
effectively three dimensional space. They become mas- 
sive as a result of the Debye screening which eventually 
determines their masses (8). These fields are dynamical 
fields, and the corresponding Maxwell term i(B^ a ^ -B^) 
is expressed in eq. (7) in terms of the cr' a - ) fields. 

We now turn to analysis of the topological density op- 
erator g(x) computed for background monopole's config- 
urations. It can be represented in terms of the effective 
scalar er field as follows 



<z(x) 



16tt 2 



rtr 
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(21) 



a=l k=l 



The vacuum expectation value of equals 2/j, a ir/NL, 
where fi a ai, = and it plays the role of the Higgs field in 
this model as explained in [12, 13]. One can explicitly see 
that the topological charge in formula (21) for a single 
monopole or antimonopole is properly normalized 3 Q = 
J d A xq(x) = ±l/N. 

Now we introduce truly singlet abelian field /jfc(x) = 



J2a=i i^-T' ) ^jk'( x ) sucn that the topological density 



3 We believe a short historical detour on fractionalization of the 
topological charge in QFT is warranted here. In given context the 
fractional topological objects appear in 2 dimensional CP N ~ X 
model [31] which were coined as instanton quarks (instanton 
partons). These quantum objects carry fractional topological 
charge Q = ±1/7V, and they are very similar to our monopoles 
in deformed QCD discussed in section II A. These objects do 
not appear individually in path integral; instead, they appear as 
configurations consisting N different 1/N objects such that the 
total topological charge of each configuration is integer. Never- 
theless, these objects are highly delocalized: they may emerge 
on opposite sides of the space time or be close to each other 
with similar probabilities. Later on, similar objects have been 
discussed in a number of papers in different context [32], [33], 
[34], [35], [36], [37]. In particular, it has been argued that the well- 
established 8/N dependence in strongly coupled QCD unambigu- 
ously implies that the relevant configurations in QCD must carry 
the fractional topological charges, see review preprint [34] and 
references therein. The weakly coupled deformed QCD model 
[12, 13] considered in this paper is a precise dynamical realiza- 
tion of this idea. 



operator g(x) for background monopoles is expressed in 
terms of this new field as follows, 
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One should emphasize that the expression on the right 
hand side of eq. (22) does not represent all properties 
of the topological density operator, e.g. it does not in- 
clude all non-abelian fluctuations which are present in 
the system. Rather it should be treated as a classi- 
cal background long distance behaviour of g(x). Fur- 
ther to this point: this new abelian field /j/c(x) is not 
an abelian projection describing the abelian U(1) N mag- 
netic monopoles in this model. Rather, the combina- 
tion e Ijfe 9j/jA;(x) describes the gauge invariant topologi- 
cal density distribution (22). Field /jfc(x) itself is not a 
gauge invariant object under generic gauge transforma- 
tions. Indeed, from relation q(x) = d^K^{x) one can in- 
fer that e yfe /jfe(x) transforms like K l (x). Transformation 
properties of K^{x) are well known: though this object 
does not carry the colour indices, it is not a gauge in- 
variant object. In fact, K^(x) transforms in a nontrivial 
way under the large gauge transformations. In particular 
/ d 3 a fl K fl determines the winding number of a "degen- 
erate" vacuum state. For our specific case of deformed 
QCD we need to know the transformation properties for 
the following operators, 

[e^Si/j^x)] transforms like q(x) (i.e. invariant) 

[e ijk f jk (x)] transforms like K\x). (23) 

Furthermore, /jfc(x) field does not discriminate different 
types of monopoles and anti-monopoles which are clas- 
sified by the affinc root of the Lie algebra a a £ A a ff. 
Instead, this new field is sensitive exclusively to the topo- 
logical charge density of these objects Q = J d A xq(yi) = 
±1/7V, not to their abelian magnetic charges ~ a a . As 
a final comment: there is no Maxwell dynamical term 
for this field. This is in drastic contrast with dynami- 
cal Maxwell term (Vc) 2 describing the abelian magnetic 
components in the effective action (7). There is no mys- 
tery here: as we shall discuss below the dynamics of 
/jfe(x) is governed by pure topological field theory with 
no Maxwell counterpart. 

One can define the gauge potential (x) in association 
with tensor field fj k (x) introduced above and new scalar 
potential a(x) describing the divergent portion of this 
tensor as follows: 

1 

2 f 

The di(x),a(x) potentials in eq.(24) are not directly re- 
lated to the abelian magnetic potentials and effective <r 
fields discussed above. The fractional topological charge 
of the monopoles can be expressed in terms of a(x) po- 
tential as follows, 



f jk (x) = [dja k (x) 



cWx)] - -e^c^x). (24) 
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where surface E defines the boundaries of our system. We 
have to take this surface to infinity if we define our system 
on R 3 . Our normalization is chosen in a such a way that 
a single monopole classified by a a and fractional topo- 
logical charge Q = 1/7V is described by an effective long 
distance field a(x) which satisfies V 2 a(x) = — 47ri5 3 (x) 
with asymptotic behaviour a(x) = l/|x|. 

Our next step is to insert the delta function into the 
path integral with field 6(x) as a Lagrange multiplier 



/ 



5 q(x) 



1 



[e«*$/ ifc (x)] 
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V[b]e^ dix 'W-tfW-raT^^'feWl), (26) 



where <j(x) ~ tr 



in this formula is treated as 



the original expression (4) for the topological density op- 
erator including the fast non-abelian gluon degrees of 
freedom, while /jfc(x) is treated as a slow- varying ex- 
ternal source describing the large distance physics for a 
given monopolc's configuration, similar to the treatment 
in ref. [17] of external currents for quasiparticles. 

Our task now is to integrate out the original non- 
abelian fast degrees of freedom and describe the large dis- 
tance physics in terms of slow varying fields in form of the 
effective action S[ar, fjk,b]. We use the same scmiclassi- 
cal approximation as we did before which is expressed in 
terms of the low energy effective action (6). The only 
new element in comparison with previous computations 
is that the fast degrees of freedom must be integrated 
out in the presence of new slow varying background fields 
fjk,b which appear in eq. (26). Fortunately, the compu- 
tations can be easily performed if one notices that the 
background field 6(x) enters eq. (26) exactly in the same 
manner as 8 parameter enters (6). Therefore, assuming 
that 6(x) is slow varying background field we arrive to 
the following effective action 

Z~ J V[b]V[(T]V[f]e- St °^ b -n- s ^' T ^ (27) 

S top[bJ] = 4^iV / 3 dSxb (x)e ijk diMx) 
d 3 x6(x)V 2 a(x); 
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2L V2tt 
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J ^ 3 t^i v 



+ 6(x) 
N 



There are two new elements in comparison with our pre- 
vious expression (7). First, the topological term S top 
emerges. This term can be also written as 



Stop = —i I d 



4tt7V 



6 4 (x) = a,6(x), (28) 



which brings 4 it into the line with conventional expres- 
sion employed in the Higgs model [17]. The second 
new element which appears in (27) is that Sduai^fr] 
now depends on pure topological field 6(x) which has 
no Maxwell counterpart. 

One should comment here that we neglected the sur- 
face terms in expression for Stop [b, /] such that only 
scalar potential a(x) from eq. (24) enters the final ex- 
pression for Stop [b, f] ■ These surface terms very often are 
crucial in similar studies in condensed matter systems de- 
fined on a finite manifold with the boundaries which may 
have nontrivial topologies such as torus. These surface 
terms are known to be responsible for the dynamics of 
the so-called "edge states" in topological field theories. 
We come back to this point later in the text. However, 
for the purpose of this work the surface terms can be ne- 
glected as we mostly discuss a trivial topology S 2 in this 
work. There is no physical degeneracy in this case and 
the system itself is characterized by a single unique vac- 
uum state. Nevertheless, the topological feature of the 
theory such as the topological long range order described 
by topological action (27), (28) does not go away when 
topologically trivial manifold is considered. 

In fact, these features are manifested in a different 
way. To be more specific, the main goal of the rest of 
the section is to argue that the well-known resolution 
of the celebrated U(1)a problem is a direct consequence 
of the topological order described by topological action 
(27), (28). In different words, we want to argue that the 
(would be) Goldstone boson receives its mass in this sys- 
tem (in apparent contradiction with conventional sym- 
metry arguments) as a result of topological features of 
the Chern-Simons action (28). 



B. Topological susceptibility in BF theory. 

We now want to compute the correlation function 
(g(x),g(0)) entering the expression for the topological 
susceptibility (15) by integrating out b, a fields 



(g(x),g(0)) = - 



1 f P[&p[<rp[a]e- s V 2 a(x),V 2 a(0) 



(AttNL)' 



To carry out the computations we limit ourself by consid- 
ering = vacuum state where (cr) = for the massive 
er fields. We expand the cos term in (27) by keeping the 
quadratic term for long range field fe(x), 



N 



o=l 



6(x) 
N 



(N 



1 - 



1 /6(x) 



N 



(29) 



4 In fact, the constraints on field 6;(x) from [17] require that 
fei(x) ~ diA when the boundary of M 3 is a topologically trivial 
S 2 , such that we are not loosing much information with identifi- 
cation 6j(x) = 9j6(x) in eq. (28). 
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Now, the obtained Gaussian integral over T>[b] can be 
explicitly executed, and we are left with the following 
Gaussian integral over V [a] 



S[a] 



1 1 



d 3 3 



2N( (4tt)2 J r3 
As the next step we rescale a(x) field 

a'(x)- aW 



to bring S [a] to more conventional form 



a(x)V 2 V 2 a(x) 



S[a'] 



(fx 



a'(x)V 2 V 2 a'(x) 



(31) 



(32) 



With this normalization, the corresponding Gaussian in- 
tegral over J T>[a'] can be easily computed 



-S[a>] 



W(x), W(0) 



-S\a> 



<5 3 (x), (33) 



where S[a'] is denned by eq. (32). Now we are ready 
to complete the computations of the topological suscep- 
tibility using the topological BF action (27). We express 
the original topological density operator (25) in terms of 
a(x) and take into account the expression for Gaussian 
integral (33). Final expression for the gauge invariant 
correlation function 



(g(x),<z(o)) 



C 



NL- 



■<S 3 (x) 



(34) 



precisely reproduces the original formula (15) which 
was derived without any mentioning of any auxiliary 
fields [23]. One can also compute a gauge variant cor- 
relation function 



lim 

k-*0 



Akx 



(Via(x),Vja(0)) 



k^ k j 



(35) 



This object is very similar to the computations (13) using 
the Veneziano ghost. The unphysical pole (35) has pre- 
cisely the same nature as the pole (13) in the Veneziano 
construction (13). In fact, the transformation properties 
of our field VjCi(x) are the same as -FQ(x) field (23) in 
the Veneziano construction (13). 

Based on this observation and comparing (13) with 
(35) we identify our topological fields constructed for the 
deformed QCD with the effective Veneziano ghost. This 
identification uncovers the nature of the Veneziano ghost 
as an effective topological non-propagating field. In both 
cases this pole is not related to any physical massless 
degrees of freedom which are not present in the system 
as the theory is gapped. Rather, it contributes to the 
non-dispersive portion of the gauge invariant correlation 
function (12), (14), (34). Still, this unphysical topological 
field does contribute to the dependent portion of the 
ground state energy. 



In weakly coupled deformed QCD one can carry out all 
the computations without even mentioning the topolog- 
ical fields or the Veneziano ghost as formula (15) shows. 
However, formulation of this phenomenon in terms of the 
topological QFT reveals its deep nature which is other- 
wise hard to understand. As explained above this contact 
term (34) is not related to any physical propagating de- 
grees of freedom. In computations (15) it emerges as a 
result of the tunnelling transitions between the degen- 
erate topological sectors. The non-dispersive nature of 
this term in present computations (based on topologi- 
cal effective Lagrangian (27)) is manifested itself in sat- 
uration of the xym by non-propagating, non-dynamical 
long range 6(x),a(x) fields. These fields are not dynam- 
ical fields as they do not have the Maxwell counterpart. 
Nevertheless, these fields are crucial as they saturate the 
non-dispersive contact term in topological susceptibility, 
as we have demonstrated above. 

Entire framework advocated in this paper is in fact a 
matter of convenience rather than necessity. The same 
comment also applies to CM systems: the BF formula- 
tion [14-19] using the topological QFT is simply a mat- 
ter of convenience to represent the known and previously 
established results (such as braiding phases, the ground 
state degeneracy, etc) using a beauty of topological quan- 
tum field theory. As we shall discuss in section IV the 
manifestations of this long range order in QCD and in 
CM systems arc somewhat different, but the beauty of 
the topological BF formulation remains the same. 



C. The mass generation for (would be) Goldstone 
boson in topologically ordered system 

We want to reproduce the behaviour (17) in deformed 
QCD using an appropriate generalization of the topolog- 
ical BF action (27) when the massless quarks are intro- 
duced into the system. This study will further illuminate 
the relation between the auxiliary a(x),6(x) fields and 
the unphysical Veneziano ghost. To proceed with this 
task we have to introduce the light matter field which 
is represented in this model by the rj -field [23]. If 
Xym were vanished the rj would be conventional massless 
Goldstone boson which is nothing but the phase of the 
chiral condensate. However, xym 7^ in 4d QCD (12) 
as well as in deformed QCD (15). In different words, the 
rj field receives its mass exclusively as a result of gener- 
ating of the nonzero contribution to Xym with "wrong 
sign" (12), which is the key element in resolution of the 
celebrated U(1)a problem [24-29]. In context of the 
present work we want to see how the rj physical con- 
tribution exactly cancels (as Ward Identities require) the 
non-dispersive term in the topological density (17). We 
want to see how it happens in BF formulation using the 
topological action (27). It will shed a new light on a 
very deep relation, already mentioned above, between the 
Veneziano ghost and the topological a(x), 6(x) fields. Es- 
sentially we want to see how the (would be) Goldstone 
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boson becomes a massive field in topological QFT in ap- 
parent contradiction with conventional symmetry argu- 
ments. 

In the dual sine- Gordon theory the rj' meson field ap- 
pears exclusively in combination with the 6 parameter as 
6^0 — <f>(x), where <j> is the phase of the chiral conden- 
sate which, up to dimensional normalization parameter, 
is identified with physical 7/ meson in QCD. As it is well 
known, this property is the direct result of the transfor- 
mation properties of the path integral measure under the 
chiral transformations ip ~^ ex p(*75^)V'- Therefore, 4>(x) 
enters the effective action exactly in the same way as the 
6(x) field does as it couples to the topological density op- 
erator exactly in the same way (26). Therefore, one can 
integrate out the fast degrees of freedom exactly in the 
same way as we did before to arrive to the following effec- 
tive action which now includes the (would be) Goldstone 
boson field 



V{b\D{(T\D[a\D[(j)]e 



-(5to P +S dual [ CT ,b,0]+S ) 



(36) 



S 4 = [ d 3 ^ (V0) 5 



S to P [&, a] = 



4ttN , 

Sdual[<X,M = / d3x ^{£) 



d d a;6(x)V 2 a(x); 
2 



N 



cos a„ ■ a 



+ b(x) - 0(x) 
N 



where coefficient c determines the normalization of the </> 
field and has dimension one. This coefficient, in princi- 
ple, can be computed in this model, but the results do 
not depend on its numerical value. In terms of these 
parameters the (f> mass is given by [23], 



<' = ON' 



(37) 



There are two new elements which appear in (36) in com- 
parison with (27). First, kinetic term for (f> field is gener- 
ated and parametrized by S<j>. Secondly, </3(x) field enters 
the .S'duail ', b, <fi] in a very precise and specific way con- 
sistent with Ward Identities. 

Our task now is to compute the topological susceptibil- 
ity in the presence of massless quark field by integrating 
out the auxiliary 6(x),a(x) fields. We repeat the same 
steps which led us to (30). The only difference is that 
a(x) and </>(x) fields arc mixing such that effective action 
SQCD[a,4>] now takes the form 



(«(x),«(o)) QCD = i J 



V[a]e- S ^ CD V 2 a(x), V 2 a(0) 



Sqcd[ci, 1 



d 3 x 



1 



2iVC(47r) 5 
|(W(X)) 2 



4irN 



a(x)V 2 V 2 a(x) 



V0(x) • Va(x) 



(38) 



where we performed the integration by parts of the mix- 
ing term J d 3 xcf>V 2 a = — J d 3 x\7cj) ■ Va. Our next step is 

to eliminate the non-diagonal term J d 3 xV(f> ■ Va in (38) 
by making a shift 



«(*) _ 



(x) + 



AttcN 



a(x), 



(39) 



where 1/y/c is inserted to the left hand side of eq. (39) 
in order to bring the kinetic term for <^(x) field to its 
canonical form. The problem is reduced to the compu- 
tations of the Gaussian integral similar to (30) with the 
only difference that the effective action now, after the 
rescaling (31), takes the form 



SqcdW \(f> 2 ] = r 



(V> 2 (x) 



d xa' (x) 



V 2 V 2 -m 2 ,V 2 



(40) 
a'(x) 



which replaces previous expression (32). In formula (40) 
parameter m 2 , is the 77' mass in this model and it is given 
by eq (37). In terms of this rescaled field a'(x) the Gaus- 
sian integral which enters (38) can be easily computed 
and it is given by 



/2?[q , ]e- s Q^V 2 a / (x),V 2 a , (0) 
J V[a']e- s QOD[a'] 



= [<5(x)-m 2 ,G(x)] (41) 



where Sqcd is defined by (40) and the massive Green's 
function G m (x) = e is normalized in conventional 



Airr 



way m 2 / d 3 xG m (x) = 1. Collecting all numerical coef- 
ficients from (31), (38) and (41) the final expression for 
the topological susceptibility in the presence of massless 
quark takes the form 



(q(x),q(0)) QCD = 



NL 2 



<5(x) 



Airr 



(42) 



It precisely reproduces our formula (17) which was de- 
rived by explicit integration over all possible monopole's 
configurations without even mentioning the topological 
auxiliary fields. It now has two terms: the first term 
with "wrong sign" which has non-dispersive nature, and 
which was computed previously (34) using the same aux- 
iliary topological fields. It also has a new dispersive 
term related to the massive 77' propagating degrees of 
freedom. It has a negative sign in accordance with gen- 
eral principles (11). The Ward Identities are satisfied 
Xqcd = J d 3 x(q(x),q(0)) Q cD = as a result of ex- 
act cancellation between the two terms. The celebrated 
U(1)a problem is resolved in this framework exclusively 
as a result of dynamics of the topological a(x) , 6(x) fields. 
These fields are not propagating degrees of freedom, but 
nevertheless they generate a crucial non-dispersive con- 
tribution with "wrong sign" which is the key element of 
formulation and resolution of the U(1)a problem and the 
generation of the 7/ mass. 
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D. Topological fields and the Veneziano ghost. 

The expression for the correlation function (41) with 
action (40) can be represented in a complementary way 
which makes the connection between the Veneziano ghost 
and topological fields much more explicit and precise. 

To proceed with out task, we use a standard trick 



V 2 V 2 



m 2 ,V 



to represent the 4-th order operator 

which enters the effective action (40) as a combination 
of two terms with the opposite signs: the ghost (f>i and a 
massive physical <j) field. To be more specific, we write 



1 



1 



1 



V 2 V 2 



m 2 ,V 2 



y 2 -m 2 v , 



such that the Green's function for a(x) field which en- 
ters the expression for the topological susceptibility (41) 
can be represented as a combination of two Green's func- 
tions: for physical massive field with conventional kinetic 
term and the ghost field with "wrong" sign for the kinetic 
term. As usual, the presence of 4-th order operator in eq. 
(40) is a signal that the ghost is present in the system. 
This signal is explicitly manifested in eq. (43). The con- 
tact term in this framework is precisely represented by 
the ghost contribution. Indeed, the relevant correlation 
function which enters the expression for the topological 
susceptibility (41) can be explicitly computed using ex- 
pression (43) for inverse operator as follows 

/ 2?[a']e- s ^[ Q 'l VV(x), VV(0) 



-Sqco[o.'] 



(44) 



d 3 p _ t 
re 



ipx 


p 4 


1 


m 2 , 




h m 2 , 


ipx 


p 2 






p2 _ 


\- m 2 , 





5{x) - m 2 , 



Airr 



which, of course, is the same final expression we had be- 
fore (41) with the only difference is that it is now ex- 
plicitly expressed as a combination of two terms: physi- 
cal massive rf contribution and unphysical contribution 
which saturates the contact term with "wrong" sign. 

A new twist in this computation is to represent the 
correlation function (44) by introducing two fields </>i(x) 
and </>(x) replacing the a'(x) which enters the effective 
action (40) as the 4-th order operator. To be more pre- 
cise, we rewrite our action (40) in terms of new fields 
0i (x) and 0(x) as follows 



S, 



QCDl 



d 3 3 



d 3 3 



V&(x)) 
(V0(x 



V<Mx) 
m 2 ,</> 2 (x) 



(45) 



while a'(x) field is expressed in terms of new fields <^i(x) 
and </>(x) as 

1 



a'(x) = — (<£(x) - 0i(x)) , 



(46) 



while the topological density q{pc) operator is expressed 
in terms of these fields as follows 



q = 



C 



NL 2 



■V a = 



NL 2 m 2 , 



(47) 



This redefinition obviously leads to our previous result 
(42), (44) when we use the Green's functions determined 
by the Lagrangian (45) for the physical massive field </> 
and the ghost 4>i, 



(q(x),q(0)) QCD = 



NL 2 



<5(x) — m 2 , - 



4irr 



(48) 



It is quite amazing that precisely this structure (45) had 
emerged previously in study of of the U{1)a problem in 
2d Schwinger model in Kogut-Susskind (KS) formulation 
[38], see also [39] with related discussions. Topological 
density operator in 2d Schwinger model e^F^ is also 

expressed as e llv F llv ~ d^d^ (<f> — (f)\j similar to (47). 

Furthermore, one can show that this structure still holds 
even with non- vanishing quark mass m q , in which case 
an additional term ~ m q cos(</> + 4>2 — 4>i) appears in ef- 
fective action (45), similar to analogous expression in KS 
description [38]. In fact, our notations for 4>,<pi,<p2 fields 
entering (45) are precisely the same notations which had 
been used in ref. [38] to emphasize the similarity. Fur- 
thermore, analogous structure also emerges in 4d QCD 
when the the topological density operator is expressed in 
terms of the Veneziano ghost where q ~ □(</> — 4>i ) has 
precisely the same structure [40] . 

An important point here is that the contact term in 
this framework is explicitly saturated by the topological 
non-propagating auxiliary fields expressed in terms of the 
ghost field </>i, similar to 2d KS ghost or 4d Veneziano 
ghost 5 . From our original formulation without any aux- 
iliary fields reviewed in section II B it is quite obvious that 
our theory is unitary and causal. When we introduce the 
auxiliary fields (which arc extremely useful when one at- 
tempts to study the long range order) the unitarity, of 
course, still holds. Formally, the unitary holds in this for- 
mulation because the ghost field 4>\ is always paired up 
with 4>2 in each and every gauge invariant matrix element 
as explained in [38] (with the only exception being the 
topological density operator (47) which requires a special 



It is important to emphasize that KS and Veneziano ghosts 
should not be confused with the conventional Fadeev-Popov 
ghost which is normally introduced into the theory to cancel 
out unphysical polarizations of the gauge fields. Instead, the 
KS, Veneziano ghost is introduced to account for the existence 
of topological sectors in the theory, see [39] for references and 
details. In four dimensional case the Veneziano ghost can not 
be confused with Fadeev-Popov ghost as the Veneziano ghost 
being a singlet does not carry a colour index, in contrast with 
Fadeev-Popov ghost. The sole purpose of the Veneziano ghost 
is to saturate the contact term with "wrong sign" in topological 
susceptibility, similar to eq. (48) in deformed QCD model. 
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treatment presented in this section). The condition that 
enforces this statement is the Gupta-Bleuler-like condi- 
tion on the physical Hilbert space H p h ys which reads like 

(02-0i) (+) |^ P hys) -0, (49) 

where the (+) stands for the positive frequency Fourier 
components of the quantized fields. The crucial point 
here is that the formulation of the theory using the topo- 
logical fields has an enormous advantage as the long range 
order is explicitly present in formulation (27), (36) and 
therefore, in equivalent formulation in terms of the ghost 
field as eq. (45) states. 

• Our arguments, based on analysis of a simplified ver- 
sion of QCD essentially suggest that this key element 
of the U(1)a problem represented by eq.(15) is a direct 
manifestation of the long distance topological action (27), 
(28), (36). Similar structure in CM systems is known to 
describe topologically ordered phases. It is naturally to 
assume that the deformed QCD also belongs to a topo- 
logically ordered phase. Furthermore, one can explicitly 
see from our computations above that the rj' generates 
its mass as a result of mixture of (would be) Goldstone 
field with topological auxiliary field. Therefore, we inter- 
pret the well known resolution of the U(1)a problem in 
deformed QCD as a result of dynamics of a topological 
field describing the long range order of the system. 

IV. PHYSICAL INTERPRETATION. 
ANALOGIES WITH CONDENSED MATTER 
SYSTEMS. 

The central line of all our previous discussions is that 
the resolution of the celebrated U(1)a problem in de- 
formed QCD is a direct manifestation of the topological 
order of the system. The main argument is based on 
analysis of the topological action (27), (28), (36) which 
exactly reproduces the crucial element of the U(1)a prob- 
lem, the topological susceptibility as cqs. (34), (42), (48) 
demonstrate. A similar topological action in condensed 
matter (CM) systems is known to lead to a variety of very 
non-trivial properties as a manifestation of topologically 
ordered phases realized in these systems, see [14-19] and 
many references therein. 



A. Differences and Similarities between CM 
systems and deformed QCD 

As we already mentioned, there is a fundamental dif- 
ference between CM systems defined in Minkowski space 
time d = (2 + 1) and Euclidean 3d "deformed QCD" 
which has been studied in the present work. In particu- 
lar, instead of propagating quasiparticles in CM systems 
we have pseudo-particles (monopoles) which saturate the 
path integral. As a result of this difference wc can not 
use many standard tools which normally would detect 
the topological order. For example, we can not compute 



the braiding phases of charges and vortices which are 
normally used in CM systems simply because our system 
does not support such kind of excitations. Further to 
this point, the topological 6(x),a(x) fields entering the 
abelian BF action (27), (28), (36) are not related in any 
way to the physical E&cM field in contrast with CM sys- 
tems where topological a M field couples to the physical 
E&cM charges. It prevents us from forming a real vor- 
tex which carries magnetic field. Also, these topological 
6(x),a(x) fields do not have canonical Maxwell counter- 
parts. Instead, the abelian topological field in our case 
decouples from E&zM charges, and behaves like field 
as discussed in section III A, see eq.(23). 

Furthermore, the "degeneracy" in deformed QCD 
model is related to degenerate of winding states \n) which 
are connected to each other by large gauge transforma- 
tion, and therefore must be identified as they correspond 
to the same physical state. It is very different from con- 
ventional term "degeneracy" in topologically ordered CM 
systems when distinct degenerate states are present in 
the system as a result of formulation of a theory on a 
topologically non-trivial manifold such as torus. 

In case of deformed QCD one should anticipate a simi- 
lar behaviour when the system is defined on topologically 
nontrivial manifold. To be more specific, if the bound- 
ary of the Euclidean space M 3 in eq.(l) is a topologi- 
cally trivial § 2 than one should expect a unique vacuum 
state. At the same time, if the Euclidean space R 3 in 
eq.(l) is additionally compactified on a large torus, i.e. 
R 3 — > T 2 x K 1 , than one should expect an additional 
topological Z2 symmetry for SU(2), in close analogy with 
behaviour of CM system [17]. In fact, the emergence of 
this additional topological Z2 symmetry has been explic- 
itly demonstrated in deformed QCD in [41]. This symme- 
try is realized as physical degeneracy in the limit of large 
T 2 . Such a behaviour of the system further supports our 
claim that the deformed QCD belongs to a topologically 
ordered phase. Indeed, a conventional gapped theory can 
not be sensitive to any variation of the boundary condi- 
tions at arbitrary large distances. 

The moral is: the systems in both cases demonstrate a 
huge sensitivity to the large distance physics formulated 
in terms of the boundary conditions. For example, even 
if we formulate a CM system on a topologically trivial 
manifold such as sphere, there would be a unique ground 
state in this "trivial" case. However, it is quite obvi- 
ous that the long range order in this system governed by 
topological BF action, does not go away. One can not use 
the conventional tools, such as degeneracy of the ground 
state, to detect a topological order. Instead, one should 
use some different observables to detect the long range 
order which obviously remains in the system even if it is 
formulated on a topologically trivial manifold. 

In fact we discussed previously some similarities be- 
tween long range order in CM systems and in QCD in 
[20]. In particular, the analysis of the Aharonov -Casher 
effect as formulated in [42] in many respects is very sim- 
ilar to our discussions in QCD. The relevant part of that 
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work can be stated as follows. If one inserts an exter- 
nal charge into superconductor when the electric field is 
exponentially suppressed ~ cxp(— r/X) with A being the 
penetration depth, a neutral magnetic fluxon will be still 
sensitive to an inserted external charge at arbitrary large 
distance. The effect is pure topological and non-local in 
nature. The crucial element why this phenomenon oc- 
curs in spite of the fact that the system is gapped is very 
similar to the key point in deformed QCD discussed in 
this section. Indeed, the system has different topological 
states u n (number of Cooper pairs) are analogous to the 
topological sectors |n) in our work. As a result of the 
"tunnelling" , an appropriate ground state U(9) must be 
constructed as discussed in [42], analogous to the \9) vac- 
uum construction in gauge theories. This state U(9) is 
an eigenstate of the so-called "modular operator" which 
commutes with the hamiltonian. In our work an analo- 
gous role plays the large gauge transformation operator T 
such that T\9) = exp(— i9)\9) which also commutes with 
the hamiltonian [T, H] = 0, such that our system must be 
transparent to topologically nontrivial pure gauge config- 
urations, similar to transparency of the superconductor 
to the "modular electric field" from rcf. [42]. 

As we already mentioned at the beginning of this sec- 
tion: there is a fundamental difference between our sys- 
tem and CM analogues because our topological fields, 
while abelian, do not have a Maxwell counterpart, do 
not couple to the physical charges, and do not propa- 
gate, in contrast with physical ESzM field discussed in 
ref. [42]. Nevertheless, both systems behave very simi- 
larly as a result of "degeneracy" mentioned above. The 
dynamics of these degenerate states can be studied in 
terms of the long range auxiliary topological field (or, 
what is the same, in terms of a la Veneziano ghost) as we 
discussed in section III. 



B. "String- net condensation" in CM [43] vs long 
range "Skeleton" in the lattice QCD [1, 2] 

While we demonstrated a number of supporting ar- 
guments that QCD indeed lies in a topologically ordered 
phase, its microscopical nature remains a mystery. In this 
section we attempt to uncover (at least partly) this mys- 
tery by presenting an analogy with CM systems where 
it has been suggested [43] that the microscopical mech- 
anism for topological phases is the so-called "String-net 
condensation" . We shall argue below that a similar struc- 
ture has been observed in the lattice QCD [1, 2] which 
was coined as "Skeleton". In fact, the lattice results 
[1, 2] was our original motivation for the present stud- 
ies in search for a long range order in QCD, sec Intro- 
duction. By providing this analogy with "String-net con- 
densation" in CM we speculate that the emergence of the 
"Skeleton" is the microscopical mechanism leading to the 
long range topological order in QCD. 

The "Skeleton" represents the extended, thin, coher- 
ent, locally low-dimensional configurations embedded in 



4d space, with opposite sign sheets interleaved. We 
present few arguments below to argue that "string -net" 
condensate in CM systems as constructed in [43] is a di- 
rect analog of the "Skeleton" structure observed in lattice 
simulations in QCD [1-10]. 

Indeed, in both cases the configurations themselves 
have lower dimensionality than the space itself. How- 
ever, these low-dimensional configurations are so dense, 
and they fluctuate so strongly that they almost fill all 
of entire space. In both cases an effective tension of 
the configurations vanishes as a result of large entropies 
of the objects which overcome the internal tension. It 
leads to the condensation of the "string-nets" and per- 
colation of the "Skeleton" correspondingly. If the effec- 
tive tensions of these configurations were not vanished, 
we would observe a finite number of fluctuating objects 
with finite size in the system instead of observed perco- 
lation of the "Skeleton" and condensation of the "string- 
nets" . Furthermore, typically a "Skeleton" is spreading 
over maximal possible distances percolating through the 
entire volume of the system similar to "string-nets" which 
condense. Finally, in both cases the V and T invariance 
holds as a result of presence of two coherent networks. 
In String -net condensation this is achieved by consider- 
ing two topological QFT's with opposite chiralities. In 
"Skeleton" studies there are two oppositely- charged sign- 
coherent connected structures (sheets). The V and T in- 
variance holds in QCD as a result of delicate cancellation 
between the opposite sign sheets. This delicate cancel- 
lation may be locally violated as a result of an external 
impact such as heavy ion collision. Apparently, such kind 
of local V and CV violation indeed has been observed in 
heavy ion collision experiments at RfflC, Brookhaven, 
and LHC, Geneva, see section VB for the details and 
references. 

The crucial difference between the two cases is of 
course the nature of the objects: CM systems live in 
real Minkowski space-time while lattice QCD measure- 
ments are done in Euclidean space-time where corre- 
sponding configurations saturate the path integral. In 
the present context it implies that while in CM sys- 
tems the corresponding "string -nets" are made of real 
particles/quasiparticles organized into extended objects 
which may condense, in QCD the corresponding ex- 
tended "skeleton" configurations live in 4d Euclidean 
space. Therefore, they should be interpreted as the ob- 
jects describing the tunnelling events in Minkowski space 
time, similar to instantons. The term "condensation" 
normally used in CM literature is not quite appropriate 
for such 4d objects. Therefore, it is more appropriate 
to describe the relevant physics by term "percolation" of 
extended configurations. 

What could be an appropriate analytical tools to study 
this structure in QCD? In deformed weakly coupled QCD 
we argued [44, 45] that the required structure may in- 
deed emerge in form of the domain walls correspond- 
ing to interpolation between topologically different but 
physically equivalent winding vacuum states. When one 
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slowly moves from weak coupling regime to strong cou- 
pling regime these extended domain walls become very 
crumpled and wrinkled objects with large number of fold- 
ings and effectively vanishing tension as a result of large 
internal entropy. In this case an arbitrary number of such 
objects can be formed and they can percolate through 
the entire volume of the system, forming a kind of a vac- 
uum condensate consisting these crumpled objects. Nev- 
ertheless, the topological charge distribution with such 
striking feature as extended coherence should still per- 
sist as the transition from weak to strong coupling regime 
should be smooth [12, 13]. It is quite likely that an ap- 
propriate description for this physics in strong coupling 
regime should be formulated in terms of holographic dual 
model rather than in terms of quantum field theory, as 
argued in [20], however we leave this subject for future 
studies. 

To conclude this section: we speculate that the "Skele- 
ton" structure observed in the QCD lattice simulations 
[1-10] may represent a basic microscopical mechanism 
supporting the topological long range order in QCD, sim- 
ilar to "string -net" condensation in CM systems [43] . It 
is clear that much work needs do be done before this spec- 
ulation becomes a workable framework. Nevertheless, in 
next section V we discuss some possible manifestations 
of the long range order in QCD if our speculations and 
analogies presented above in this section IV turn out to 
be correct. 



V. POSSIBLE MANIFESTATIONS OF THE 
LONG RANGE ORDER IN QCD 

In section III we argued that the resolution of the cel- 
ebrated U(1)a problem is a direct manifestation of the 
topological order of the system. The arguments were 
based on study of the weakly coupled deformed QCD 
model. In section IV we argued that the long range struc- 
ture observed in the lattice simulations [1-10] may be 
considered as an evidence that the transition from weakly 
coupled deformed QCD to the strong coupled QCD is in- 
deed smooth as similar "strings- net" structure in CM 
is normally attributed to topologically ordered phase. If 
the topological order indeed persists in this transition 
to the strong coupling regime than we should expect a 
number of other profound consequences (along with res- 
olution of the celebrated U(1)a problem which has been 
already discussed) of this long range topological order in 
QCD.' 

The main goal of this section is to argue that we might 
be indeed witnessing some of the manifestations of this 
long range order. We elaborate below on two phenom- 
ena which are well established experimentally. At the 
same time, a theoretical understanding of these effects 
is still lacking. We present some arguments suggesting 
that these two phenomena can be naturally understood 
if QCD indeed belongs to a topologically ordered phase. 



A. Universal hadronization temperature 

T H ~ (160 - 170) MeV 

Over the years, hadron production studies in a variety 
of high energy collision experiments from e + e~ annihila- 
tion to pp and pp interactions with energies from a few 
GeV up to the TeV range, the production pattern always 
shows striking thermal aspects with apparently quite uni- 
versal temperature around Th — (160 — 170) MeV. It is 
very difficult to understand the nature of this "apparent 
thcrmalization" as one can not even speak about kinetic 
equilibration. In different words, the thermal spectrum 
in all high energy collisions emerges in spite of the fact 
that the statistical thermalization could never be reached 
in those systems. Hagedorn concluded that the hadrons 
are simply born in equilibrium [46], in apparent contra- 
diction with causality. 

We argue below that the source for such striking ther- 
mal features could be related to the coherent structure 
of sheets making the "skeleton", see section IV B. In this 
case the "apparent thermalization" emerges as a result 
of tunnelling events (described by these coherent 4d con- 
figurations) accompanied by particle production, rather 
than a result of interaction of the produced particles. 

As we discussed in section IV B the "skeleton" con- 
figurations should be interpreted as tunnelling processes 
which arc happening all the time in Minkowski vacuum 
with no interruptions. These tunnelling processes in 
Minkowski vacuum (when no any external sources are 
present in the system) do not lead to any emission or 
absorption of real particles, similar to the persistent tun- 
nelling events in Bloch's case. These tunnelling events 
simply select an appropriate ground state of the system 
which is a specific superposition of winding |n) states. 

However, when an external source emerges in the sys- 
tem (e.g. as a result of high energy collision) a delicate 
cancellation mentioned in section IV B between the op- 
posite sign sheets will lead to a production of real parti- 
cles. We expect that the high energy collisions including 
e + e~, pp and pp interactions do not completely destroy 
the coherent vacuum structure discussed in section IV B 
as the vacuum remains in the same confined topologi- 
cally ordered phase. However, the collisions lead to some 
distortions of these low-dimensional sheets rather than 
breaking them. It is crucial that these Euclidean long 
range vacuum configurations with vanishing width (as 
observed in lattice simulations) actually describe instan- 
taneous tunnelling events (with no violation of causality) 
in Minkowski space. Relaxation of these long range con- 
figurations to their conventional form (which existed in 
vacuum before the collision occurred) leads to emission of 
real physical particles. 

What is expected spectrum of the emitted particles as 
a result of this mechanism? With a number of additional 
assumptions one can explicitly see emergence of the ther- 
mal spectrum as a result of the Unruh-Hawking radia- 
tion [47]. We refer to papers [40, 48, 49] and references 
therein where all relevant formulae in the present con- 
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text with emerging Planck spectrum have been discussed. 
One should emphasize that the Planck spectrum in this 
approach is not resulted from the kinetics. Rather, the 
Planck spectrum in high energy collisions in this frame- 
work is resulted from the stochastic tunnelling processes 
when the event horizon emerges as a result of strong in- 
teractions. In such circumstances the spectrum must be 
thermal with T given by 

T=^-~k QC D- (50) 

The acceleration a in this framework describes the dy- 
namics of distorted low-dimensional sheets which orga- 
nize the "skeleton" as discussed in [49]. This parameter 
in case of a small distortion (such as peripheral heavy ion 
collisions, see section VB) can be, in principle, expressed 
in terms of the auxiliary axion field, similar to a(x), 6(x) 
fields from section III, see footnote 7 to avoid confusion 
with term "axion" . In general, computation of a is a 
hard problem of strongly coupled QCD. 

An approximate universality of the temperature with 
no dependence on energy of colliding particles nor their 
nature (including e + e~, pp and pp collisions) in this 
framework is due to the fact that the emission occurs 
from the distorted long range vacuum configurations 
rather than from the colliding particles themselves. Fur- 
thermore, it is expected that the acceleration a (as a 
dynamical field in this framework) is not fluctuating on 
AqJ, d scales, but sensitive to larger scales such as size 
of the system. This is because its dynamics is governed 
by long range topological fields entering the topological 
action (similar to (27), (28) in deformed QCD) rather 
than by conventional fast gluon fluctuations. In deformed 
QCD this property can be easily seen from the action 
written in form (45) when topological fields arc repre- 
sented by the ghost <\>\ (x) which does not "know" about 
Aqcd scale. It should be contrasted with physical mas- 
sive field 0(x) which of course "knows" about Aqcd as 
to,,' ~ Aqcd ■ The initial energy of colliding particles de- 
termines the normalization of the fluctuating topological 
field, but not its spectral properties. 

What could be an appropriate technique to study this 
question in strongly coupled regime? As suggested in 
[20, 48], it is very likely that the dual holographic de- 
scription might be an appropriate technical tool to an- 
alyze these questions. However, we leave this subject 
for future studies. Much work needs to be done before 
this proposal initiated in [48] and further elaborated in 
[40, 49] becomes a workable framework in strongly cou- 
pled regime. 

However, we believe that the crucial element of this 
idea (which is formulated in terms of the tunnelling 
events of the long range coherent "skeleton" described 
by the auxiliary axion field, similar to a(x),6(x) fields 
from section III) has a potential to eventually explain 
Hagedorn's notion that the states were prepared before 
the collision occurs. This explanation would not violate 
causality as topological auxiliary fields are not propagat- 
ing degrees of freedom in this framework. In different 



words, these fields can not send or receive signals as they 
describe entanglements rather than correlations. 



B. Local V and CP violation in QCD. Chiral 
Magnetic Effect (CME) 

The violation of local V and CV invariancc in QCD 
has been a subject of intense studies for the last couple of 
years as a result of very interesting ongoing experiments 
at RHIC (Relativistic Heavy Ion Collider) [50, 51] and, 
more recently, at the LHC (Large Hadron Collider) [52- 
55]. 

The main idea to explain the observed asymmetries is 
to assume that an effective 8(x,t)i n d ^ is induced as 
a result of collision. As the 0(x,t)i n d parameter enters 
the Lagrangian (6) with the topological density operator 
q(x), the local V and CV invariance of QCD is broken on 
the same scales L where 9(x, t)i n d ^ is correlated. As a 
result of this violation, one should expect a number of V 
and CV violating effects (such as Chiral Magnetic Effect 
(CME), Charge Separation Effect(CSE), Chiral Vortical 
Effect (CVE)) taking place in large region L as suggested 
in [56]. 

It is important to emphasize that 9(x, t)i n d must be 
correlated on large scales L ^> Aq^ D ; otherwise, all V 
and CV odd effects will be washed out. Furthermore, 
the effective Lagrangian approach advocated in [56] is 
justified if and only if 9{x,t) in d is correlated on large 
scales L ^> Aq CD . These ideas with specific applications 
to heavy ion collisions were further developed in [57, 58] , 
see recent review papers [59, 60] with a complete list of 
references and latest developments in this active field, 
including some recent experimental works. 

One of the crucial questions for the applications of the 
CME to heavy ion collisions is a correlation length of the 
induced (9(x,t)i n d) ^ 0. Why are these V odd domains 
large? In applications to heavy ion collisions the size of 
the system L ~ 10 fm is much larger than conventional 
Aq^ ~ 1 fm. Apparently, a relatively large correlation 
length is a required feature for interpretation of the ob- 
served asymmetry [50-55] in terms of CME as argued in 
[40, 45, 49]. 

We think that the crucial element in understanding 
this key question is deeply rooted to the properties of 
topological order as we advocate in this work. To be 
more concrete: as we discussed above, the emission of 
real particles in this framework is resulted from relax- 
ation of the distorted long range sheets of the "skeleton" 
to its vacuum configuration. These distortions (resulted 
from the collision) spoil the delicate exact cancellations 
between oppositely charged long range vacuum configu- 
rations which existed in vacuum before the collision oc- 
curred. The size of these V odd distortions is order of 
the size of the system, i.e. ~ L. This long range or- 
der explains why CME, CSE and CVE are operational 
in this system and how the asymmetry can be coherently 
accumulated from entire system. Exactly this long range 
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order explains how (9(x, t)i n d) 7^ could be correlated on 
large scales. It also justifies the assumption made in [56] 
where 0(x,t)i n d was treated as an external background 
field in effective Lagrangian approach with correlation 
length much larger than any conventional QCD fluctua- 
tions, L 3> Aq CD . 

The key element in estimation of measured asymme- 
tries is algebraic sensitivity to acceleration parameter a 
introduced above 6 . This parameter for peripheral col- 
lisions with large impact parameter can be estimated as 
a ~ L _1 where L _1 is the size of the system. The al- 
gebraic sensitivity to large distances leads to sufficiently 
large intensity for asymmetries as estimated in [49]: 

asymmetries from [50-55] ~ —7 — - — ~ 10 -3 , (51) 

LAq CD 

where electric charge e ~ y/a ~ 10 _1 enters the estimate 
(51) because CME, CSE and CVE effects are propor- 
tional to anomalous coupling of the topological fluctua- 
tions with EkM field. Also, the factor « ~ iV -2 - 0.1 
in eq. (51) describes conventional suppression ~ iV~ 2 
related to any 9 dependent effects. Finally, for numerical 
estimates we assume Aqcd/s ~ LAqcd — 10. If a long 
range order were not present in the system, we would get 
a strong suppression cxp(— AqcdL) instead of power like 
suppression (LAqcd) 1 hi eq. (51). 

Numerical estimation (51) is consistent with intensi- 
ties of the observed asymmetries presented in refs. [50- 
55]. Furthermore, this entire framework is consistent 
with all known experimental results. In particular, the 
correlations due to the local V violation should demon- 
strate the universal behaviour similar to the "univer- 
sal apparent thermalization" discussed in section V A as 
the source for the both effects is the same and related 
to long range "skeletons" as argued above. For exam- 
ple, the effect should not depend on energy of collid- 
ing ions. Such independence on energy is indeed sup- 
ported by observations where correlations measured in 
Au+Au and Cu+Cu collisions at ^/snn = 62 GeV and 
\J snn = 200 GeV are almost identical and independent 
on energy. The same tendency also continues for the LHC 
energies: the ALICE Pb-Pb results for 3-particle correla- 
tor at ^/satjv = 2.76 TeV are almost identically coincide 
with RfflC Au+Au results at y/sJfN = 200 GeV, which 
is precisely what was anticipated as the sizes for Au 79 
and Pb 82 are almost the same, and therefore, parameter 



The algebraic sensitivity with respect to a can be tested using 
the Veneziano ghost model[40, 61, 62]. As we discussed in sec- 
tion III the topological fields a(x), b(x) in deformed QCD essen- 
tially represent the ghost field governed by the effective action 
(45). Separately: the structure (45) exactly coincides with corre- 
sponding expression in 4d QCD when U(1)a problem is treated 
using the Veneziano ghost. Therefore, the testing of the algebraic 
sensitivity using the Veneziano ghost in [40, 61, 62] essentially 
implies that the same algebraic sensitivity to a will emerge in de- 
formed QCD when the topological fields a(x), 6(x) are expressed 
in terms of the ghost field. 



L entering (51) is also the same, see [49] and references 
to the experimental results therein. 

The algebraic decay (a/AQco) p with p ~ 1 which ap- 
pears in estimate (51) is a direct manifestation of the 
topological long range order in QCD. This algebraic scal- 
ing can be, in principle, experimentally tested as the size 
of the system entering (51) is related to the size of the 
colliding heavy ions. Some suppression of the measured 
correlations with increasing the size of the system indeed 
has been observed. For example, the effect for Au+Au 
collisions with A ~ 197 is suppressed in comparison with 
Cu+Cu collisions with A ~ 64, see [49] and references 
to the experimental results therein. 

We conclude this section with the following proposal 
based on analogy discussed in section IV. We argued that 
the long range order in CM systems and in QCD is re- 
sulted from very similar structures: "strings -net" con- 
densation in CM and "skeleton" in QCD correspondingly. 
We also argued that in both cases the V and T invariance 
holds on average as a result of presence of two coher- 
ent networks with delicate cancellations between them. 
In "string -net" condensation the delicate cancellation 
is achieved by considering two topological QFT's with 
opposite chiralities. In "skeleton" studies there are two 
oppositely- charged sign-coherent connected structures 
such that all V and T odd effects vanish on average. 
As we discussed above in heavy ion collisions (in event 
by event studies) this cancellation is not perfect and it is 
not absolute, and in principle may lead to some V odd 
effects as a result of formation of large V odd domains 
which can be described in terms of the auxiliary axion 
field 7 . This field is similar to our topological 6(x) field 
in our studies of the deformed QCD in section III, see 
cqs. (27), (28). 

It would be interesting to find a similar effect in CM 
systems where a local violation of V and T invariance 
may also occur as a result of non-cancellation between 
two topological QFT's with opposite chiralities as a result 
of some external impact which locally may destroy its del- 
icate cancellation. The point is that this non-cancellation 
may happen on a macroscopical scale of a sample, rather 
than on a scale of the correlation length determined by a 
gap of the system. In fact, the effect similar to CME and 
CSE[56] (mentioned above in regards to heavy ion colli- 
sions) has been discussed in CM literature where it was 
coined as topological magneto-electric effect [64]. Fur- 
thermore, the topological effects considered in [64] are 
described in terms of the domain walls which are very 
similar (in spirit) to the domain walls studied in weakly 
coupled deformed QCD in [44, 45], and which were inter- 
preted there as a direct manifestation of the long range 



7 We use term "axion" which is a jargon here. There is no real 
new dynamical degree of freedom such as axion in QCD, see 
reviews [63] about physical axion as real degree of freedom and 
dark matter candidate. However, the V odd domain described in 
terms of 9i n d{x) justifies our terminology as the dynamical 0(x) 
is, by definition, the axion. 
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structure of the system described by a coherent configu- 
ration of the axion field (6i n d) ■ 



VI. CONCLUSION AND FUTURE 
DIRECTIONS 

The main "technical" result of the present work can be 
formulated as follows. The analysis of a simplified version 
of QCD suggests that the non-dispersive contribution to 
the topological susceptibility (12) (which is a key ele- 
ment of formulation and resolution of celebrated U(1)a 
problem) emerges because the deformed QCD can be de- 
scribed in terms of the auxiliary topological a(x),6(x) 
fields governed by BF-like action (27), (28), (36). Such 
an effective action in CM systems is normally consid- 
ered as a signal for emergence of a topologically ordered 
phase with a number of known striking features. There- 
fore, it is naturally to assume that the deformed QCD 
also lies in a topologically ordered phase. While we can 
not use many standard tools (such as the braiding phases, 
physical degeneracy etc as explained in section IV) which 
would confirm the topological order, we still can feel the 
long range order in the system through different phe- 
nomena discussed in this paper. One could argue that 
this long range order persists in strongly coupled QCD 
as well, as there should be no any phase transitions in 
passage from "deformed QCD model" to real QCD. 

However, the computations in "deformed QCD model" 
do not say much about the source, the nature, the "mech- 
anism" of this long range order. At this point we return 
to our first line in introductory section I on very puzzling 
recent lattice results [1-9] which actually motivated this 
study. What are the observable manifestations of the 
long range order observed in these simulations? One can 
inverse the question: what is the microscopical realiza- 
tion of the topological order in QCD if it is confirmed 
by future studies? In CM systems it has been shown 
that a physical mechanism for topological phases can be 
formulated in terms of the string- net condensation [43]. 
We argued in section IV B that the structure observed 
on the lattices and coined as "skeleton" is analogous to 
"String-net condensation" in CM systems. In different 
words, we argued that the microscopical mechanism for 
the long range order in QCD if it is confirmed by future 
numerical and analytical computations could be precisely 
the long range structure which already had been studied 
[1-9]. 

However, the observable manifestations of this topo- 
logical ordered phase is much more difficult to detect as 
explained in IV B because EhM field as a probe does 
not couple directly to the topological structure observed 
on the lattice simulations in contrast with CM systems 
where external Eh,M field is a perfect probe of the long 
range order. Nevertheless, there are different manifesta- 
tions of the long range order such as generating the mass 
of would be Goldstone bosom (the so-called U(1)a prob- 
lem in QCD) as well as some other manifestations which 



were shortly reviewed in section V. We think that all 
these phenomena are already signalling us that the long 
range order is present in the system. Much work needs 
to be done before the speculations reviewed in section V 
become a working framework. 

A first step along this line of study could be testing the 
ideas formulated in section VB using some CM systems, 
which arc known to belong to a topologically ordered 
phase. Specifically, a local violation of V and T invari- 
ance (similar to observed at RHIC and the LHC) may 
occur in a CM system formulated in terms of "string-net 
condensation" of two topological QFT's with opposite 
chiralities. Such a local violation could be a result of 
non-cancellation between two topological QFT's due to 
some external impact which may locally destroy the del- 
icate cancellation between the two "string-networks" . It 
would be interesting to see if this non-cancellation may 
occur on a macroscopical scale of a sample, rather than 
on a scale of the correlation length determined by a gap 
of the system. It would be also interesting to test in CM 
systems the idea that "would be Goldstone" boson may 
generate its mass through mixing with a topological field, 
similar to what happens in deformed QCD as discussed 
in section III C. 

Finally, one should note that the crucial physics which 
is responsible for a number of striking features discussed 
in this paper and which are due to the tunnelling events 
between "degenerate" winding states 8 can be in principle 
tested in the Maxwell theory defined on a compact man- 
ifold [65]. To be more specific, when the Maxwell system 
is defined on a compact manifold there will be tunnelling 
events which are very similar in nature to the effects dis- 
cussed in the present work. Corresponding topological 
effects are also not related to the physical propagating 
degrees of freedom, similar to our discussions of the topo- 
logical susceptibility in section III. Furthermore, the re- 
sulting topological contribution to the Casimir pressure, 
while numerically very small has the sign which is oppo- 
site to the well known expression for the Casimir effect, 
similar to the properties of the non-dispersive contribu- 
tion to the topological susceptibility discussed in section 
III. Essentially, this extra topological contribution is pre- 
cisely the source of a tiny violation of a commonly ac- 
cepted (but generally wrong) receipt that the Casimir 
effect due to Maxwell photons could be obtained by mul- 
tiplying the corresponding scalar expressions by a factor 
of two. This extra contribution is exactly analogous to 
the topological terms discussed in this work, which can 
not be expressed in terms of conventional propagating 
degrees of freedom. 



see a comment on term "degeneracy" in the given context in 
footnote 2. 
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